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ABSTRACT 

We use the accretion model to explain the early exit of rj Car from its 2009 X-ray minimum. In the 
accretion model the secondary star accretes mass from the primary wind near periastron passage, a process 
that suppresses the secondary wind. As the shocked secondary wind is responsible for most of the X-ray 
emission, the accretion process accounts for the X-ray minimum. The early exit from the 2009 X-ray minimum 
after four weeks, instead of ten weeks as in the two previous minima, is attributed to the primary wind that 
during the last minimum was somewhat faster and of lower mass loss rate than during the two previous X-ray 
minima. This results in a much lower mass accretion rate during the X-ray minimum. We show that using 
fluctuations in these quantities that are within the range deduced from fluctuations in the X-ray flux outside 
the minimum, can account for the short duration of the last X-ray minimum. The shorter X-ray minimum 
may have further implications on the recovery of the system from the spectroscopic event. 

Subject headings: (stars:) binaries: general— stars: mass loss— stars: winds, outflows— stars: individual {rj Car) 



1. INTRODUCTION 

For 2023 days we have been waiting for the intrigu- 
ing periastron passage of the massive binary system rj 
Car. And indeed, we were not disappointed, as the 
new X-ray observations (Corcoran 2009) brought some 
surprises. The periodicity itself, which is observed in 
all wavelengths (e.g., radio, Duncan & White 2003; IR, 
Whitelock et al. 2004; visible, van Genderen et al. 
2006, Fernandez Lajus et al. 2008; UV, Smith et al. 
2004; emission and absorption lines, Nielsen et al. 2007 
Damineli et al. 1997, 2008a, b; X-ray, Corcoran 2005, 
2008, Hamaguchi et al. 2007), is attributed to a binary 
orbital period (Damineli 1996). 

The X-ray light curve was followed by RXTE (Corco- 
ran 2009) during three minima. In all three minima the 
X-ray intensity increases prior to the start of the X-ray 
minimum and then drops sharply to minimum, which 
was about the same level in the three minima (Corcoran 
2005, 2009). However, while the 1998 and 2003.5 min- 
ima lasted 65 — 70 days (Corcoran 2005), the exit from 
the last X-ray minimum started only 27 days after peri- 
astron passage (Corcoran 2009); last periastron passage 
is taken at January 11, 2009. 

In Kashi & Soker (2009a; see also 2009b) we have al- 
ready pointed out the possibility that the duration of 



the X-ray minimum can be shorter or longer. We con- 
sidered relatively small variations in the primary wind 
properties, and estimated that the X-ray minimum can 
be different by several days. In this Letter we reexamine 
X-ray variations far from the minimum, i.e., near apas- 
tron, and deduce that the variations (fluctuations) in the 
primary wind velocity and mass loss rate can be up to 
a factor of two We then show (fj3]) that within the 
accretion model (Soker 2005; Akashi et al. 2006; Kashi & 
Soker 2008, 2009a) such an increase in the primary wind 
velocity near periastron passage can account for the early 
exit of 77 Car from the X-ray minimum. In the accretion 
model the X-ray minimum is attributed to the accretion 
of the primary wind by the secondary star. The accretion 
process suppresses the secondary wind, that otherwise is 
responsible for most of the X-ray emission (Corcoran et 
al. 2001; Pittard & Corcoran 2002; Akashi et al. 2006; 
Henley et al. 2008; Okazaki et al. 2008; Parkin et al. 
2009). Our discussion and summary is in ^ 

The present study presents neither new ingredients to 
the accretion model nor new type of calculations. The 
essence of this paper is to strengthen the accretion model 
by showing that it can easily accommodate the new ob- 
servations by only changing the variables that are most 
likely to change from orbit to orbit and along the orbit, 
namely, the primary wind properties. 
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2. POSSIBLE FLUCTUATIONS IN THE PRI- 
MARY WIND 

The binary parameters are as in our previous pa- 
per (Kashi & Soker 2009a, where references are given): 
The assumed stellar masses are Mi = I2OM0 and 
AI2 — 3OM0, the eccentricity is e = 0.9, and the or- 
bital period is P = 2024 days. The stellar winds' 
mass loss rates and terminal velocities that were used 
are Afi = 3 x IO-^Mq yr'^, M2 = lO'^M© yr-\ 
vi^oo = 500 km s-i and W2,oo = 3000 km s'^. While 
in our previous papers the mass loss rate and terminal 
velocity of the primary wind were held constant with 
the above values, in the present paper we vary them to 
account for the early exit from the X-ray minimum. 

To estimate the largest fluctuations in the primary 
wind properties, we use the fluctuations in the X-ray 
luminosity near apastron in the following manner. The 
hard X-ray emission observed in 1] Car is emitted by the 
shocked secondary wind (Corcoran et al. 2001; Pittard 
& Corcoran 2002; Akashi et al. 2006; Henley et al. 2008; 
Okazaki et al. 2008; Parkin et al. 2009). For constant 
secondary wind properties the X-ray luminosity varies as 
~ , where D2 is the distance of the stagnation point 
of the colliding winds from the secondary star. This can 
be understood from two different considerations. When 
the radiative cooling time is long, the X-ray emission is 
(Akashi et al. 2006) ~ 0.5M2swi(T/2/rcooi2), where 
Tf2 — D2/V2 is the outflow time of the shocked gas, and 
^00012 oc n^^ is the radiative cooling time of the shocked 
secondary wind. The postshock secondary wind density 
varies as n cx D^^. We find therefore oc D^^. In the 
second approach the X-ray emissivity (power per unit 
volume) An^ is used. The luminosity is = V2An^, 
where V2 oc is the volume of the shocked secondary 
wind. The mass in the volume V2 is proportional to its 
outflow time r/2 — D2/V2. Therefore, n oc -D^^ (as 
before), and again we recover the relation Lx oc D^^ . 

The distance D2 is given by equating the momentum 
fluxes of the two winds. When 77 Car is near apastron 
(hence orbital velocity is negligible) D2 is approximately 
given by D2 — ?'C/(1 + 0? where here we define C, = 
{M2V2/ MiViY^"^ ^ and r is the orbital separation. For 
our typical parameters C — 0.45. Using this relation in 
the expression for the X-ray luminosity gives 

Lx = K{r,V2,M2)^^, (1) 

where the function K depends on the orbital separation 
and secondary wind properties. The logarithmic varia- 
tion of the X-ray luminosity with respect to the primary 
wind momentum discharge pi = MiVi is given by 

dLx_ _ 1 dpi 
Lx " 2(1 + Pi ■ 



From the X-ray light curve (Corcoran 2005) we find 
that the variations between cycles and during a short 
time within one cycle can be dL^/Lx — ±0.25. If we 
attribute this variation to changes in the primary wind 
properties, then we find from equation ([2|) and C = 0.45 
that dpi/pi ~ ±0.7. Basically, if the mass loss rate does 
not change the terminal velocity can vary in the range 
^ 300 — 850 km s~^; such changes are seen in the solar 
wind. If the wind velocity is larger when the mass loss 
rate is lower, then the variation in the wind velocity can 
be even larger. In the acceleration zone close to the pri- 
mary star the variations can be much larger. In such an 
eccentric, asynchronous binary system like rj Car, tidal 
interaction between the stars can result in kinetic en- 
ergy dissipation through viscous shear. This can lead 
to variations in mass-loss rate and wind velocity struc- 
ture, or even asymmetric mass-shedding (Koenigsberger 
& Moreno, 2009). 

Our conclusion is that variations in the primary wind 
velocity by a factor of < 2 in the acceleration region are 
reasonable. As mentioned earlier, in our previous paper 
(Kashi & Soker 2009a) we assumed variations in pi of 
no more than ~ 10%, and predicted the possibility of an 
early (or late) exit by several days. Motivated by the 
new observations (Corcoran 2009), we reexamined the 
primary wind variations, and concluded in this section 
that we can allow for much larger variations. 

3. VARYING THE ACCRETION RATE AT 
THE LATEST PERIASTRON PASSAGE 

In Kashi & Soker (2009a) we took two extreme pro- 
cesses, that we expect to bound the true mass and an- 
gular momentum accretion rates, to estimate the mass 
and angular momentum accretion rate. These are the 
Bondi-Hoyle-Lyttleton (BHL) accretion process from a 
wind, and a Roche lobe overflow (RLOF) type mass 
transfer. Very close to periastron passage, \t\ < 10 day, 
the accretion process is an hybrid of the BHL and the 
RLOF mass transfer processes, but at the end of the 
accretion phase the accretion will be of the BHL type. 
The calculated accreted mass in the different models 
considered was estimated to be in the range Mace ~ 
0.4 — 3.3 X 10~^Mo, and the accretion rate was typically 
Mace = 5 X 10"'^ - 5 X 1O"^M0 yr"\ with the higher 
values close to periastron. 

The accretion processes close to periastron is very 
complicated, and for its accurate study one must use 3D 
hydrodynamical numerical codes. In this highly eccentric 
binary system the primary angular velocity is expected 
to be much below the angular velocity of the secondary 
near periastron. Therefore, the mass transfer process will 
not be as in the RLOF process in synchronized binary 
systems. For that reason, in the present paper we use 



only the BHL accretion process. 

We take the /3-profile to describe the primary wind 
acceleration 

/ R \^ 

vi{ri) ^ Vs + {vi.oo - Vs) \1 - y-j , (3) 

where Ri is the primary radius, ri distance from the 
primary center, and Vs = 20 km s~^. For our "standard 
case" we use the parameters /3 = 3 and a terminal wind 
velocity of wi,oo — 500 km s"-'^. The orbital eccentricity 
is taken to be e = 0.9. 

The mass accretion rate for the standard case was 
calculated in our previous paper (Kashi & Soker 2009a). 
The BHL mass accretion rate as calculated there is drawn 
by the solid very thick blue line in Figure [TJ At t — 
65 days, where the system starts leaving the 1998 and 
2003.5 X-ray minima, the mass accretion rate is 

Mref = M^ccit = 65 days, /3 = 3, vi,oo = 500 km s^^) = 

8.7 X IO^'^Mq yr"^ 
(4) 

This is depicted by the dashed horizontal thin blue line 
in Figured] We take this value to be the reference value 
for the accretion rate at which the secondary wind re- 
builds itself after the accretion phase. We note that for 
the present goal the exact value of Mj-of is not impor- 
tant. What is important is how we vary the primary 
wind properties to obtain the same value at a much early 
time. 

We emphasize that in our model for an early exit rel- 
ative to the previous two cycles, it is required that the 
primary wind velocity in the last minimum was higher 
than the wind velocity in the previous two minima. It is 
not required that the primary wind velocity after perias- 
tron be higher than the primary wind velocity far from 
periastron. In the present paper we take the velocity 
of the primary wind in the previous two minima to be 
vi = 500 km s^^. 

We now search for the value vi,ooNow in the /? = 
3 profile that would give an accretion rate of 8.7 x 
10~^Mo yr~^ at i = 27 days, the time when the early 
exit from the 2009 X-ray minimum started. We keep the 
mass loss rate at AIi — 3 x lO^^M© yr~-^. We find that 
fi.ooNcw — 900 km s^^ gives the desired accretion rate, 
as shown by the solid thick red line in Figure [TJ This is 
within the reasonable range of the primary wind velocity 
fluctuations we estimated in the previous section. 

We also find that if we keep vi^oo = 500 km s^^ and 
Ml = 3 X 1O~^M0 yr~^ as in the standard case, then 
a value of /3 = 0.35 gives the desired accretion rate at 
< = 27 days. This case is presented by the solid thin 
green line in Figure [1] 




t [days] 



Fig. 1. — The variation of accretion rate near periastron passage. 
Only the Bondi-Hoyle-Lyttleton (BHL) accretion prescription is 
presented. In the first ~ 20 days extra accretion is possible via a 
RLOF type mass transfer. The different lines, from top to bottom, 
depict the mass accretion rate for the following cases characterized 
by primary mass loss rate in units of Mq yr"'^, the value of /3 in the 
acceleration profile given by equation and the primary wind 
terminal velocity give in units of km s~^ , (Mi , /3, vi^ao)'- Solid very 
thick blue line: The standard case that was used in our previous 
paper (3 X 10"*, 3, 500); Solid thick red line: (3 X 10""', 3, 900); 
Dashed thick red Une: (1.5 X 10"*, 3, 780); Dashed thin green line: 
(1.5 X 10-"*, 1.05, 500); Solid thin green line: (3 X 10"'', 0.35, 500). 
The dashed thin blue line is the reference accretion rate as given 
by equation Bt . 

It is likely that a more efficient acceleration process 
occurs when the mass loss rate is lower. We assume 
now that the primary mass loss rate was lower in the 
2009 minimum, which lead to a more efficient acceler- 
ation. To demonstrate the feasibility of our model, we 
take the mass loss rate to have been half its typical value 
Ml. 

New — 1-5 X 10 ^Mq yr ^. For Mi^ncw and keeping 
/3 = 3, we find that wi,ooNow2 = 780 km s~^ gives the 
reference accretion rate Mref at t = 27 days (dashed 
thick red line in Figure [J). For Afi^Ncw and keeping 
vi,oo = 500 km s~-^, we find that /3 = 1.05 gives the ref- 
erence accretion rate at t = 27 days (dashed thin green 
line in Figure [1]). 

As well, our model might account for the change in 
the rate of the exit from the two events. We quantify 
the observed exit rate by defining the average value of 
AL^/ At over the exit from minimum time period. Here 
At is the exit time interval between the last epoch in 
which the X-ray luminosity has been at its minimum 
value and the first epoch in which the X-ray luminosity 
has returned to its quiescent level, and ALx is the dif- 
ference in the corresponding X-ray luminosities. Using 
the RXTE results (Corcoran 2005; 2009), wc find that 
the (averaged) observed exit rate of the 2003.5 minimum 
was AL^/At = 1.92 x 10^* ergs"^. The observed exit 
rate of the 1998 was lower by 11%, while the observed 



Cycle 


2003.5 


2009 


Observations 


ti [days after periastron] 


64.7 


26.3 


t2 [days after periastron] 


91.0 


64.7 


L^i [ erg s-i] 


1.9 X 10^4 


1.9 X 10^4 


1^2 [ erg s"^] 


6.3 X 10^4 


6.5 X 10^4 


AL^/At [ erg s'^] 


1.92 X 1028 


1.39 X 1028 


Ratio (2003.5/2009) 


0.72 


Model 


(Macc)i [Mq yr-1] 


10-6 


10-6 


(Macc)2 [Mo yr-1] 








AMacc/Ai [Mq yr-2] 


-1.4 X 10^5 


-9.5 X 10^6 


Ratio (2003.5/2009) 


0.68 



Table 1: Comparing the observed exit rates ratio to the one ob- 
tained from our model, if we assume that the X-ray exit rate is 
proportional to the change in the mass accretion rate (although in 
reality no simple relation is expected). 



exit rate of the 2009 minimum was substantially lower, 
with AL^/At = 1.39 x lO^* ergs'^, which is 72% of the 
observed exit rate of the 2003.5 minimum. 

We check whether the change in time evolution of the 
mass accretion rate during the exist from the X-ray mini- 
mum might explain the change in the observed exist rate. 
For that, we measured the average rate of change in the 
mass accretion rate, AMacc/Ai, as given in Figured] 
The average is over the exit time interval of each event. 
We summarize the results in Table [T] We find that the 
ratio in the average values (AMacc/Ai) during the 2009 
and the 2003.5 minima is 0.68. The proximity of the ratio 
0.68 (ratio of modelled accretion rate in the two minima) 
to 0.72 (ratio of observed X-ray flux time evolution in the 
two minima) should not be given a strong weight, as we 
do not expect a simple relation between mass accretion 
rate and X-ray luminosity. What is important is the gen- 
eral qualitative change. We conclude that our model can 
account, with the same change in parameters, to both 
the 2009 early exit from the minimum, and for the 2009 
slower rate of increase in the X-ray flux during the exit 
from minimum. 

4. DISCUSSION AND SUMMARY 

Our aim was to explain the early exit of rj Car from 
its X-ray minimum in the last cycle. Whereas in the pre- 
vious 2 cycles the exit occurred ~ 10 weeks after perias- 
tron passage, in the last cycle it occurred only 4 weeks 
after periastron passage. It was known that in the ac- 
cretion model variations in the primary wind mass loss 
rate and/or velocity can cause variations in the duration 
of the X-ray minimum (Kashi & Soker 2009a). However, 
the exit was earlier than anticipated in our earlier paper. 
The reason is that in our previous paper we considered 



relatively small variations in the primary wind proper- 
ties. 

In the accretion model the secondary star accretes 
mass from the primary wind near periastron. While for 
most of the ~ 5.54 yr orbital period the secondary grav- 
ity has a negligible effect on the colliding winds cone, 
near periastron the shocked primary wind is very close 
to the secondary star and very dense, and the secondary 
gravitational field becomes non-negligible. As was shown 
in previous papers (Soker 2005; Akashi et al. 2006; Kashi 
& Soker 2008a, 2009a) the gravitational field of the sec- 
ondary ensures accretion. The accretion process is as- 
sumed to shut down the secondary wind, hence the main 
X-ray source; the system then enters the X-ray minimum. 

Although rj Car had an early exit from the X-ray min- 
imum in the last cycle, the beginning of the minimum 
was as in the previous two cycles, and occurred shortly 
before phase zero. The onset of the accretion phase re- 
quires an almost extinction of the secondary wind, and 
therefore a high mass accretion rate that occurs just be- 
fore periastron. During the phase interval 0.01 to 
—0.01 (—20 < t <^ 20 days) the two stars are very close 
(r < 3 AU), and the mass accretion rate is an hybrid of 
a RLOF and an accretion from a wind, with the RLOF 
type most likely the dominant process very close to peri- 
astron passage (Kashi & Soker 2009a; after t ~ 20 days 
the accretion from the wind is the relevant process, and 
for that it was used in the paper). Therefore, the pri- 
mary wind velocity does not influence much the mass 
accretion rate during the onset of the X-ray minimum. 

The above explanation holds for small and moderate 
fluctuations in the primary wind properties. The very 
early exit from the X-ray minimum in the last cycle sug- 
gest that the changes in the primary wind properties 
were large. Therefore, for our model to work we require 
than the changes occurred near periastron passage, most 
likely due to the strong tidal interaction. Most likely, the 
tidal interaction amplifled a small internal change in the 
wind properties, e.g., as might results from magnetic ac- 
tivity. Koenigsberger & Moreno (2009) calculations show 
indeed that tidal interaction can cause large changes in 
the wind properties of LBV stars in binary systems. 

Following recent RXTE observations by Corcoran 
(2009), we reexamined the behavior of the X-ray emis- 
sion when the system is far from the X-ray minimum, 
i.e., it is near apastron. Based on the behavior of the 
X-ray emission, in section [2] we found indeed that the 
variations (fluctuations) in the mass loss rate and/or ve- 
locity of the primary wind can be by up to a factor of 
2. 

Armed with this finding, we examined the possible 
implications for the accretion model. We assumed that 
the secondary wind, and hence the X-ray emission, is 



recovered when the accretion rate is as in t = 65 days 
(after periastron) in the previous cycles. The case with 
the parameters used to explain the previous 2 cycles is 
termed "the standard case" . We found that by increasing 
the primary wind velocity by a factor < 2, and possibly 
with lowering the primary wind mass loss rate, we can 
have the mass accretion rate at t = 27 days as in t = 
65 days in the previous cycles. We examined several 
cases, as drawn in Figure [1] 

It should be emphasized that we did not add any new 
ingredients to the accretion model, nor did we change 
the type of calculations. The new addition is simply 
repeating the calculation from our previous papers with 
different primary wind properties. As discussed, these 
properties are likely to vary by a large factor from orbit 
to orbit and along the orbit. This new part shows that 
the accretion model is quite robust. 

One implication of our finding is that in some lines 
that come from the primary wind during the minimum 
the radial velocity in the last minimum was larger than in 
the previous minima. We emphasize that what matters 
to the early exit from the minimum is the primary wind 
velocity in the last cycle compared with the previous two 
minima, and that the wind velocity in the minimum can 
be lower than the average wind velocity far from the min- 
imum. We note also that for accretion the velocity of the 
wind expelled in and near the equatorial plane matters. 
The wind toward our line of sight can have somewhat 
different value. We observe the binary system at inclina- 
tion angle of 41° (Smith 2006), therefore such variations 
might not be as large, as tidal effects are stronger in the 
equatorial plane." 

The shorter X-ray minimum might have further impli- 
cations. As the secondary accreted less mass than in the 
previous two cycles, the accretion disc (or belt) formed 
around it may be thinner, and its dissipation time would 
be shorter. Therefore the recovery time of different spec- 
tral lines and electromagnetic bands from the minimum 
can change. 

In the same way as the last X-ray minimum was 
shorter, future X-ray minima might be longer than 
70 days. The point is that any fluctuation in the wind, 
even a moderate one, can change the duration of the 
X-ray minimum. 
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